Lifshitz transitions and crystallization of fully-polarised dipolar Fermions in an 

anisotropic 2D lattice 
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We consider a two dimensional model of non-interacting chains of spinless fermions weakly coupled 
via a small inter-chain hopping and a repulsive inter-chain interaction. The phase diagram of this 
model has a surprising feature: an abrupt change in the Fermi surface as the interaction is increased. 
We study in detail this meta-nematic transition, and show that the well-known 2~ — order Lifshitz 
transition is the critical endpoint of this first order quantum phase transition. Furthermore, in the 
vicinity of the endpoint, the order parameter has a non-perturbative BCS-like form. We also study 
a competing crystallization transition in this model, and derive the full phase diagram. This physics 
can be demonstrated experimentally in dipolar ultra-cold atomic or molecular gases. In the presence 
of a harmonic trap, it manifests itself as a sharp jump in the density profile. 

PACS numbers: 03.75.Ss, 71.10.Hf, 71.10.Fd, 71.27. +a 



I. INTRODUCTION 

The study of phenomena which deform the Fermi sur- 
face in both electron and cold atom systems has gained 
much popularity recently. For example, it has been 
sugested that the Fermi surface shape-changing Pomer- 
anchuk instability^— may describe experiments in heavy 
fermions^— quantum hall devices,— and ruthenates , 10 ! 11 
leading to a plethora of theoretical papers on the 
subject ^22, 

The Pomeranchuk instability, which breaks rotational 
symmetry but no translational symmetries, is essentially 
a transition to an electronic nematic phased— An ana- 
logue in spin systems has also been studied^ Incorpo- 
rating other electronic analogues of liquid crystal phases 
into this picture has been put forward as a general pic- 
ture of strong correlations 2 ^ 3 , with evidence for smec- 
tic phases being observed experimentally in manganites 24 
and cuprates^ 

These electronic liquid crystal phases also have a 
strong relation to dimensional crossover phenomena, 
where one can ask the question whether an array of 
one-dimensional chains (Luttinger liquids) coupled by 
a weak inter-chain hopping t± remains strictly one- 
dimensional (confinement)^ or becomes a quasi-one- 
dimensional Fermi liquid (deconfinement). Indeed, cal- 
culation methods such as self-consistent perturbation 
theory 27 -^ and functional renormalization 29 support the 
idea that the warped Fermi-surface is unstable for suffi- 
ciently small t ± , in principle therefore leading to a Fermi- 
surface modifying transition at some finite value of inter- 
chain hopping. While the issue of the Luttinger liquid 
to Fermi liquid crossover/transition is not yet fully re- 
solved, it is possible to ask a much simpler question: 



what happens if an array of one-dimensional Fermi liq- 
uids are coupled together? A rather specific example in 
this direction was the study of coupled edge states (chi- 
ral one-dimensional Fermi liquids) in superlattices which 
exhibit integer quantum Hall effect^ 

While such a toy model may not be realistic for any 
real materials, advances in laser trapping and cooling 
technology have led to the rapidly expanding field of 
trapped ultra-cold atoms, which in the context of con- 
densed matter physics can be thought of as a sort of 
quantum analogue simulation of a bulk systemj 31 i with 
unprecedented control over disorder and interactions. It 
was suggested a few years ago^ - — that exploiting the 
dipole interaction between cold polar molecules or highly 
dipolar atoms allows further control over effective in- 
teractions, in order to build exotic strongly correlated 
phases. Of particular interest in the present context is 
the case when the atoms or molecules are fermions. 33 
On the experimental front, there has been much recent 
experimental progress towards this goal using highly po- 
lar 40 K 87 Rb molecules^ and a fermionic isotope of the 
highly magnetic atom 163 Dy. 37 On the theoretical front 
there has been a flurry of activity.— ~— 

Unlike the long- wavelength scattering induced by Fes- 
hbach resonances, dipolar interactions have a power-law 
dependence on the distance between the interacting par- 
ticles and a non-trivial dependence on the relative posi- 
tion of the two particles and orientation of the magnetic 
dipoles. In the presence of a strong polarising field the 
latter translates into a strongly anisotropic interaction 
which leads to a spontaneous (though not symmetry- 
breaking) deformation of the Fermi surface.— Indeed de- 
pending on the strength of the dipolar interaction addi- 
tional, symmetry-breaking (Pomeranchuk) Fermi surface 



2 



deformations may also occur 4i Even more interestingly, 
such polarised dipolar gases can in theory be combined 
with optical lattices to generate non-trivial tailor-made 
effective Hamiltonians. 

In this connection the present authors showed^ that 
a quasi-one dimensional (quasi-lD) optical lattice could 
be used to create a system whose phase diagram features 
Fermi liquid, stripe and checkerboard ground states, as 
well as a meta-nematic quantum phase transition into a 
state with distorted Fermi surface. The model of Ref. 
features chains within which there are no interactions, 
achieved by the alignment of the polarising field at the 
'magic angle' to the tube direction. For this particular 
orientation of the field, the interactions between particles 
on different chains are purely repulsive. An experimental 
realisation of this model would thus furnish an example 
of the coupled one-dimensional Fermi liquids mentioned 
above. Interestingly, the spontaneous Fermi surface dis- 
tortion encountered in this model corresponds also to an 
interaction-induced change of dimensionality, from quasi- 
lD (open Fermi surface) to fully 2D (closed Fermi sur- 
face) behavior. An important property of this model is 
that these various transitions can all happen in different, 
well separated, regions of phase space. Thus, each of the 
phase transitions can be studied independently without 
its properties being masked by the other ones. 

More recently other authors have studied a closely- 
related model, featuring continuum tubes rather than 
discrete chains 4£ In this case only the strictly ID limit 
was considered but on the other hand the polarising field 
was allowed to point in any direction. This leads to a 
rich variety of effective interactions and a correspond- 
ingly rich phase diagram: in addition to density wave and 
meta-nematic phases, different superfluid ground states 
are expected. 

The present work extends the theories in Refs. I39ll40l 
by considering finite temperature phase transitions in 
addition to the quantum phase transitions discussed to 
date. We also develop a fully analytic theory of the meta- 
nematic transition in the neighborhood of its critical end 
point. 

The Hamiltonian of the model that we analyze in detail 
consists of spinless fermions hopping along parallel chains 
(labelled by n), with a weak hopping and an interaction 
between nearest neighbor chains: 



i,n 

- *-L (4,n c i,n+l + cj.n+l^.") + V Pi,nPi,n+lj (1) 



where the density 



Pi.n — c in C 



(2) 



The plan for the rest of the paper is as follows: In 
Section II we will study the meta-nematic transition in 
full detail, explaining its origins and its link to the Lif- 
shitz transition. In Section III we will look at the other 



competing instability of the model to a crystalline phase. 
The complete physical picture is presented in Section IV. 
We finish with a discussion of both the specific proper- 
ties of this model, and how these may be generalized to 
a more realistic microscopic model in condensed matter 
systems. 



II. THE META-NEMATIC TRANSITION 

We first consider a phase transition closely related to 
the Lifshitz transition^ 



A. Self-Energy in self-consistent Hartree-Fock 
approximation 

We study how the interactions renormalize the spec- 
trum, within the diagrammatic perturbation theory ap- 
proach. To do this, we will introduce the Matsubara 
Green function 

G(k,iw n )=/ dTe^ T ((T T c k (T)cl)), (3) 



where T T is the (imaginary) time ordering operator, Ck("r) 
is the annihilation operator in the Heisenberg representa- 
tion evolving in imaginary time, and ((...)) is the thermal 
average. In the above expression and in what follows we 
have set Boltzmann's constant ks = 1 and Planck's con- 
stant divided by 2w, h = 1. Within the diagrammatic 
perturbation theory, the Green function can be built up 
from the non-interacting Green functions 
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lUJ r , 



(4) 



Here, 



£o(k) = — 2£ii cos k x — 2t± cos k y — ji (5) 



is the bare dispersion. The full Green function is then 
given by a Dyson series to be 



G(k, iw n ) 



1 



iuj n - f (k) - £(k,wj„) 



(6) 



where E(k, iuj n ) is the self-energy. 

We will calculate the self-energy to lowest order, which 
corresponds to the Hartree and the Fock terms: 




The first diagram is the Hartree term. It has no de- 
pendence on k x , k y or iu> n and can be trivially evaluated 
to give 



-'H 



-VN/n 



(7) 
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where N is the numer of particles and £1 is the volume 
of the system. This term simply renomalizes the Fermi- 
energy and will be ignored for now. 

The second (Fock) term has more structure. The dia- 
gram corresponds to the expression 

£ F (k,u;)=T£ |^ Go(cMen)F (k-q) (8) 

where we have introduced the Fourier transform of the 
interaction 

V(p) =2V cos ( Py ). (9) 
Now, the Matsubara summation can be performed to give 

T E G °(q,-n) = 1 + e L)/T ="(q) ( 10 ) 

where n is the Fermi function (with respect to the non- 
interacting part of the Hamiltonian) . We see first that 
TiF is independent of energy, and is given by 

E F (k x ,k y ) = 2V J ^ n(q x , q y ) cos(k y - q y ). (11) 

Note that T,p(k x , k y ) = T,p(k y ) so the Hartree-Fock self- 
energy does not depend on k x . This is a direct conse- 
quence of the interaction potential being k x -independent: 
V(k x , k y ) = V(k y ) and is a specific feature of our partic- 
ular model. In the absence of any term breaking the 
symmetry k y — > —k y (which would give rise to sponta- 
neous currents) we have n(q x ,q y ) = n(q x , ~q y ), so 

£f(M = 2V cos k y -^-r JL n(q x ,q y )cos(q y ). (12) 
J Ztt Zn 

Now, from Eq. [51 we see that the role of an energy in- 
dependent self-energy is simply to renormalize the bare 
energy spectrum: 



£(k) = &(k) + E(k). 



(13) 



Specifically, the effect of this diagram is to renormalize 
the transverse hopping to 



t*> =t ± + V 



n(qx,q v )cos(q y ). (14) 



To make this Hartree-Fock theory self-consistent, we need 
to take n(q x , q y ) to be the Fermi function with the renor- 
malized spectrum 



n(q) = 



1 



1 + e «(q)/T 



(15) 



where 



While the self-consistent expression for the renormal- 
ized transverse hopping, Eq. [TJ] looks innocent enough, 
its solution has a number of surprising features for the 
anisotropic tight binding model spectrum in question. 
Our strategy for the rest of this section proceeds in two 
steps. First we will analyze the solutions of Eq.[14]at zero 
temperature, where we will discover that far from evolv- 
ing smoothly, the renormalized hopping undergoes a first 
order jump at some critical V. We refer to this Fermi- 
surface shape changing transition as the meta-nematic 
transition, and we will examine its properties showing 
that it is the finite interaction version of the well known 
Lifshitz transition, which as we shall see is its quantum 
critical end point. 

It is in order here to clarify that the term meta-nematic 
transition is used in direct analogy with the term meta- 
magnetic transition, extensively used in recent years i 47 i 48 
In the first case the transition is caused by the exter- 
nal tuning of the bare tunneling strength whereas in the 
second case the transition (or crossover) is caused by 
the tuning of an external magnetic field. Thus in both 
cases we are dealing with the non symmetry-breaking 
version of a symmetry-breaking instability (the Pomer- 
anchuk instability in one case and fcrromagnetism in the 
other) in the form of a first-order transition tuned by a 
symmetry-breaking field (namely, lattice anisotropy and 
an externally- applied magnetic field, respectively). 

In the second step, we will proceed to finite tempera- 
ture, where we will see that this line of first order tran- 
sitions ends in a critical end-point at a certain T c . We 
will also show that this transition is a density of states 
related effect, and show why we expect it to be robust 
against higher orders in perturbation theory. 



B. The meta-nematic transition at zero 
temperature 

In Fig. [T] we plot a numerical solution of Eq. [TJ] It is 
clearly seen that there is a first order jump in the renor- 
malized hopping t* at some value of V . In order to find 
the location of this jump, the meta-nematic transition, 
we must find which solution is stable in the regime where 
there are three possible solutions. This can be done by 
computing and analyzing the total free energy of the sys- 
tem. 

The non-interacting and interacting parts of the free 
energy are 



<*i>) 



d 2 k 



&(k)n(k) 



(17) 



d 2 k 



T I {n(k) In n(k) + (1 - n(k)) ln(l - n(k))} , 

2i, r A'i, 



£(q) = —2*|| cos q x - 2*1 cos q y - fi. (16) (Fi) = - / 4^ / 4^7^(q)"(k) [1 - nCk + q)] , 

2 J (2ny J {2-kY 



For a variational derivation of these self-consistency 
equations see Ref. |39j. 



(2tt) 2 J (2tt) 2 

with the total free energy obtained by adding them to- 
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FIG. 1: [Color online] Left panel: The dashed curve shows 
the saddle point of Eq. 1141 for a typical value of fi = and 
t±/t\\ — 0.1. There is an unstable region where there are three 
solutions of this equation. To determine the true ground state 
solution, we take cuts through this and look at the energy 
(right panel) for values of V/t\\ = 4.2,4.4,4.6,4.8,5.0 top to 
bottom. The local minima correspond to the saddle point 
solution, while the true minima allow us to draw on the left 
the true renormalization of t* , (solid line in first panel) which 
has a first order jump at some critical V. 

get her 

(F) = (F ) + (F 1 ). (18) 

It is worthwhile noting that the equation for the renor- 
malization of the transverse hopping Eq. [TJ] is exactly 
the solution of Jp- = 0. We defer a discussion of higher 
order terms in the interaction to section V. 

In the right panel of Fig. [T] we plot how the total en- 
ergy changes with t* when every other parameter remains 
fixed. For V below the transition, there is a single energy 
minimum. Above the lower transition line, a local mini- 
mum at higher t* appears (the middle line in the 3- valued 
parts of Figs. [T]is always a local maximum) . Between the 
two transition lines, this local minimum becomes lower 
than the previous minimum - and there is a true first 
order jump in the value of t* . Above the upper critical 
line, the local minimum at lower t* vanishes. 

We are now in a position to calculate numerically the 
energy at its minima and find the lowest one. This al- 
lows us to determine and show in Fig. [5] where the actual 
transition takes place. In addition, we plot the size of 
the jump. Note that the jump diminishes as the value 
of V where it occurs tends to 0, with a quantum critical 
end-point at V = 0. 

While V is the natural parameter to vary in order to 
understand the nature of the transition, much insight can 
be gained by plotting the renormalized t* as a function of 
the chemical potential \x for fixed bare interchain hopping 
t± and interaction strength V (in addition, in a cold atom 
realisation it will be easier to tune t± than V). This is 
done in Fig. [3] Note the very pronounced dependence of 
t* ± on n near the transition, even for V = tn, in spite of 




FIG. 2: [Color online] The renormalized hopping as a function 
of the interaction strength for n — 0.0 (left) and /j, — — 1.5iii 
(right), and various values of the bare interchain hopping t±, 
which can be read off the value of t* ± when V = 0. The lower 
panels also show the size of the jump At*^ as a function of 
the position V at which it takes place. 




o.i - 



I , 1 , 1 , 1 , 1 

-2 -1.5 -I /, -0.5 

FIG. 3: [Color online] The variation of the renormalized hop- 
ping t* with at a fixed bare hopping t± = 0.1t|| . The differ- 
ent lines are for different interaction strengths, from bottom 
to top, V = 0, 1, 2, 3, 4 in units of . 

the first-order jump itself being very small for such small 
value of the interaction strength. 

We can now plot the actual zero-temperature phase di- 
agram, showing the single line of first order phase tran- 
sitions. This is shown in Fig. HI The 2 and 1/2-order 
Lifshitz transition^ can be regarded as the quantum crit- 
ical end-point of our meta-nematic first order transition, 
which we will discuss in more detail in the next subsec- 
tion. 



C. A theory of the quantum phase transition 

So far, all of the results have been obtained by nu- 
merical evaluation of the integrals. We now go one step 
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FIG. 4: [Color online] The upper panel shows the critical 
interaction strength V c as a function of the bare interchain 
hopping t± for four different values of the chemical poten- 
tial (from top to bottom fi/t\\ = 0.0, -0.5, -1.0, -1.5). The 
dashed lines in each of these cases bound the region where 
the energy landscape has more than one local minimum. The 
lower panel shows the size of the jump At* ± as the lines in the 
upper panel are crossed. 



further, and look at an analytic solution around the quan- 
tum critical end point (the Lifshitz transition), where the 
critical V c — > 0. This occurs at the point where the (bare) 
Fermi surface just touches the edge of the Broullin zone, 
t_L = t? PT = t\\ -j«/2. We write x = {f L -t < j FT )/t\\ and 

= (t± ~~t± T )/t\\j noting that in this notation, x will 
always be negative. By a direct expansion for small x of 
the energy, Eq. 1181 at zero temperature, we find (we have 
relegated the details of the derivation to the appendix.) 



E otx z 



In Id - 



1 



— 2x [x — xq — aV] In | 



-Vb: 



; ln 2 |d 



(19) 

where a and b are (positive) constants, depending on the 
chemical potential /x, and we have ignored any terms in- 
dependent of x. The logarithms reflect the presence of 
logarithmic singularities in the density of states (van- 
Hove singularities) at the Fermi level when x = 0. 40 
We note that a similar expansion emphasizing the non- 
analyticities around the quantum critical end point has 
been made in a related two dimensional modelj 4 ^ where 
interactions also drive the Lifshitz transition first order. 
The Energy function, Eq. 1191 can now be minimized with 
respect to x: 



dE 
dx 



oc 2 [In |x| + 1] [-{x -x ) + V(a- bx In |x|)] . (20) 



The turning points |^ = at x = ±l/e (when the first 
bracket in Eq. [20] becomes zero) are unphysical and an 
artifact of our small x expansion Eq. 1191 The remaining 
local minima of [TO] occur when 



It is worthwhile noting that when V = 0, the only solu- 
tion of this equation is x — xo, as it must be, and in fact 
it is this condition that strongly constrained the form of 
the energy functional, Eq. [TO] Now, at finite V, we know 
that the meta-nematic transition occurs when the energy 
of the two local minima becomes equal, so that the value 
of x which minimizes E jumps from one minima to the 
other. This occurs at the point xq + aV = 0, or 



V c = -Xo/a 



(22) 



where the energy functional Eq.[TOIbecomes an even func- 
tion of x. This linear relation between the critical V and 
the distance from the QPT agrees prefectly with the nu- 
merical results in Fig. [4] Now, substituting the value of 
V c into Eq.l2"Tl gives us the solution x = ± exp {—l/bV c }, 
or 



M* ± = 2e~^c 



(23) 



xq = V [a — bx In |a;| 



(21) 



This non-perturbative expression for the jump in the 
renormalized transverse hopping, which acts as an ef- 
fective 'order parameter', again agrees well with the nu- 
merics in Figs. [5] and @] We note that the expression 
is of the same form as the well-known one giving the fi- 
nite binding energy of a Cooper pair in the presence of 
and arbitrarily weak attractive interaction. Its meaning 
is that the 2 and 1 /2-order Lifshitz transition turns first 
order under the effect of interactions in a non-analytic 
way, making it an idealisation never truly realised in any 
real (and therefore, interacting) system. Note that re- 
cent experimental evidence in favour of a Lifshitz phase 
transition in Na^CoC^ seems to suggest that the tran- 
sition is indeed of first order in this material.— Other 
interaction-induced non symmetry-breaking phase tran- 
sitions involving a change of t opology o f the Fermi surface 
have been discussed in Refs. 12.49.5lll52l. 



D. Finite temperature 

In this section we study the effects of temperature on 
the meta-nematic transition. The equation for the renor- 
malization of t± is the same as for the zero temperature 
case, but with the zero temperature distribution replaced 
by the finite temperature Fermi distribution. In this case 
we perform the integrations numerically. A representa- 
tive cross section of the results are plotted in Fig. [5] As 
expected, by increasing temperature, the size of the first 
order jump gets smaller, until it finally vanishes at a 2nd 
order critical point at some (T e , V c ). 

We find the position of the jump within this unsta- 
ble region by free energy considerations identical to the 
zero-temperature case. This is plotted in Fig. [5] as a 
function of temperature T for different values of \i and 
the bare hopping t±. We note that unlike the quantum 
case, the classical second-order phase transition (where 
the jump vanishes at finite T) has a regular free energy 
expansion as the van-Hove logarithms are smoothed over 
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FIG. 5: [Color online] The solution of the saddle point 
Eg. 1141 at finite temperature (from outside to inside T/t\\ = 
0.0,0.1,0.2,0.3,0.4,0.5) for a typical value of /x = and 
t±/tu — 0.05. As temperature is raised, the size of the re- 
gion with more than one local minima decreases, and eventu- 
ally disappears altogether when we reach T c , the second-order 
critical end-point of the meta-nematic transition. 




FIG. 6: [Color online] The critical line in the V — T 
plane, for an example value of chemical potential n = 
0. From right to left, the bare hoppings are t±/tu = 
0.0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55. The 
dashed lines enclose the region of interaction strengths V for 
which the free energy has more than one local minimum. In- 
set: The size of the jump At^/tii as we cross the critical line, 
as a function of temperature T. From outside to inside, the 
values of t x /t\\ = 0.0, 0.1, 0.2. 



by temperature - leading to the standard analytic mean- 
field behavior whith the expected threshold temperature 
dependence: 



At* ± oc (T c - T) 



1/2 



(24) 



III. THE DENSITY WAVE TRANSITION 

The model, Eq. Q] also has an instability to a density 
wave (DW) state, which has broken translational sym- 
metry characterized by a non-uniform density: 



Pdw{x, y) = po + Acos(2k F x + ny). 



(25) 



This arises from the quasi-one dimensional nature of the 
model, which has an almost nested Fermi surface (which 
becomes perfectly nested in the limit t± — > or /.i = 0). 
While the physics of this transition is well established, it 
is important to investigate it in more detail in this case, in 
order to determine the regions in parameter space where 
the new meta-nematic transition is not eclipsed by this 
DW transition. 

We will estimate the transition line between the DW 
state and the 'normal' phase via the Random Phase 
Approximation (RPA). We note that this apporach is 
valid only for determining where the DW instability takes 
place, and not to describe the ordered state beyond the 
instability line. For that, one could add a translational 
symmetry breaking mean field and redo the Hartree-Fock 
theory. This would have the advantage of furnishing a de- 
scription of the symmetry broken phase as well. However, 
the presence of a DW incommensurate with the lattice 
creates many calculational difficulties, and the instabil- 
ity condition obtained in this way is equivalent to that 
obtained within the RPA. Furthermore, the use of RPA 
emphasizes the two-particle nature of the DW instability, 
in contrast to the one-particle nature of the MN transi- 
tion; thus each are controlled by different properties of 
the Fermi surface (see section V). 

We start with the density-density correlation functions 
in the Matsubara time representation: 



X(k,r) = i((T T p(k,r) iO t(k ) 0))). 



(26) 



In the non-interacting case, this is given by the 
particle-hole bubble: 

X°(k,w) 

- -/(^/s GO<k+q ' e+ " )G0(q ' e) - (27) 

Treating the interactions via the RPA then gives us: 



X(k, w ) 



X°(k,co) 



l + V(k y )X°Qz,oj)' 



(28) 



by 



The momentum dependent interaction V(k y ) is given 



V(k y ) = 2V cos k y . 



(29) 



In particular, setting k v = ±7t, this is negative, and could 
lead to an instability. In fact, the Stoner criterion is 
exactly looking for poles in response functions at zero 
frequency, u> — 0. Looking at Eq. [551 we see that this 
happens when 



which is shown as the inset to Fig. [5] 



2V"X°(k,w = 0) = 1. 



(30) 



The imaginary part of Xq vanishes at w = 0. Taking the 
real part gives the Lindhard formula for density-density 
response: 



3?X (k,o;) 



<iq n(q) — n(k + q) 
(27r)2 w -^(k + q)+e(q)' 



(31) 



where n(k) is the Fermi- function. We can include the 
effects of the Hartree-Fock self-energy in this expression 
by taking £(k) to be the renormalized energy spectrum, 
given by Eqs. HH and [T1J1 rather than the bare tight- 
binding spectrum. The most important point is that 
so long as the renormalized inter-chain hopping remains 
small, X (k, 0) is only very weakly dependent on k y . 

In general, the Lindhard function (and hence the 
Stoner criterion for the phase boundary) must be eval- 
uated numerically, however there are some limits where 
we can gain insight analytically. Let us consider first the 
case where the renormalized interchain hopping t* re- 
mains zero. Therefore, near the Fermi surface, the nest- 
ing condition £(k+Q) « — £(k) is satisfied for the nesting 
vector Q = (2&f,7t). This gives 



A (Q,w = 0) 



P(e-F)hi 



(32) 



where p(c_f) is the density of states at the Fermi surface, 
and €f ~ 2t" + fx is the chemical potential as measured 
from the bottom of the band. Thus the Lindhard func- 
tion has a logarithmic divergence at wavevector Q at zero 
frequency and zero temperature. This means that the 
ground state is indeed a density wave for any finite in- 
teraction strength, with a critical temperature estimated 
from the Stoner criterion: 



-l/p(e F )V 



(33) 



This means that with zero interchain hopping, the 
metanematic transition of the previous section will never 
be seen, as it will always be preempted by an ordering 
transition. 

For finite (but small) interchain hopping, the per- 
fect nesting is lifted (so long as we stay away from the 
particle-hole symmetric point /i — 0) and there is no 
longer a singularity: 



WtQ.^ = o.r = 0) - ,,u, \\u ( pi 



meaning that a finite interaction 



1 



pMhr(ff) 



(34) 



(35) 



must be obtained before the instability to the density- 
wave state occurs. This gives us a region of parameter 
space where as the interaction is turned up, the metane- 
matic jump takes place before the DW transition, so the 
metanematic transition can be seen. 





FIG. 7: The first and second order vertex corrections 



To go beyond these simple limits requires a numerical 
evaluation of the Lindhard function. In the next sec- 
tion, we will plot the phase boundary of the DW tran- 
sition alongside the phase boundary of the metanematic 
transition allowing the region in phase space where the 
metanematic transition 'wins' to be clearly seen. Before 
going on to this however, we will take a look at correc- 
tions to the RPA, which in our particular model will be 
proven to be small. We will then look at a strong coupling 
expansion which will complement the weak coupling per- 
turbative analysis in this section, before finally giving a 
completely different derivation of the DW instability for 
the special case of t± = 0. 



A. Vertex Corrections 

In this subsection we go beyond the RPA and study the 
vertex corrections (note the fact that RPA becomes ex- 
act for a strictly ID linear spectrum^). First we examine 
the case where the renormalized hopping remains zero, 
so the Green functions are independent of transverse mo- 
mentum, and the interaction lines are independent of the 
longitudinal momentum. 

The first order correction is shown in Fig. [7] (a). We 
note that the particle lines are dressed by the Hartree- 
Fock self-energy, see Eq. 7. This takes care of all first- 
order vertex corrections except for that represented by 
the diagram (a), which is given by 



rfq 



dp 



de 2 
2^ 



Xl ^ k y^' J ( 27r )2 J 2n J (2tt) 2 
xG°(k x +q x ,ei + uj)G°(k x + q x + p x ,ei + e 2 + w) 
xG (q x ,e 1 )G (q x +p x ,e 1 +e 2 )V{p y ) = 0. (36) 

The result is zero due to the integration over transverse 
momentum, J dp y V(p y ) — 0. We emphasize that this is 
a special feature of the model under consideration. 

The second order diagrams are also shown in Fig. [7] 
Diagrams (b)— (e) will all be zero for the same reason as 
the first order diagram. Diagram (/) is therefore the first 



8 



correction to the RPA 

K x, y, ) J (2?r)2 J 2n J (2?r)2 J 2jt 

G°(k x + q x ,e 1 + u))G°(k x + q x + p x ,ti + e 2 + 
G°(q x , e 1 )G°(q x + p x ,e, + e 2 )X°( Pyi e 2 )[V ( Py )] 2 (37) 

While this is non-zero, note firstly that the entire ex- 
pression is independent of k y , it is smaller than X° by 
a factor (V/t^) 2 and most importantly, does not have 
a divergence at k x — 2kp (basically the divergence is 
cured by the extra integrations). In the region of interest 
k x w 2kp, lj = 0, it can therefore be considered as a small 
perturbation to X° within the RPA and thus does not 
have to be taken into account to get the essential physics 
correct. 

All of this is strictly true for t± = 0. Nevertheless, it is 
easy to check that when t± ^ 0, these vertex corrections 
form a power series in t±/tn. Therefore as long as the 
renormalized value of t± remains small, the RPA is a 
very good approximation. Similar features has been seen 
e.g. in a model of stacked edge states in the case of the 
integer quantum Hall effect. 30 

B. Strong coupling limit 

While the smallness of the vertex corrections means 
that RPA is formally a very good approximation for the 
model in question, it still remains a weak coupling per- 
turbation expansion. We can complement these results 
by showing that in the strong coupling limit, the ground 
state is indeed a (2fcj?,7r) density wave. This demon- 
strates the qualitative correctness of the phase diagram 
we will present in the next section, even if the RPA does 
not estimate with high precision the position of the tran- 
sition line for large interaction strengths. 

We can proceed by first taking the kinetic part of the 
Hamiltonian Eq. [TJto zero t\\ = t± = 0. In this limit, we 
can rewrite the Hamiltonian as: 

Ho = V ^ Pi,nPi,n+l (38) 
i,n 

= \ V ^2 [(Pi,n + Pi,n+l) ~ (Pi,n - Pi,n+l) 2 ] 
i.n 

where we have used the property pf n = pi >n . From the 
above form of the Hamiltonian it becomes obvious that 
a uniform distribution of fermions over the chains is not 
energetically favorable. To illustrate this, consider the 
half-filled case, where half of the sites are occupied in the 
system with M chains each one having N c sites, therefore 
total fermions N = MN c /2. If the real space pattern is 
such that every column (perpendicular to chains) has ev- 
ery second site occupied then the total energy is zero. 
Similarly this result stays the same if we choose any par- 
ticular way how to "combine", or "lock" these columns 



and form the chains. Therefore there is a macroscopic 
degeneracy of 2^"=. 

It is interesting to map the problem onto an Ising 
Hamiltonian by allowing the values of a pseudospin S 
to be +1 in the presence of a fermion at a given lattice 
site and -1 in the absence of it. Then we can transform 
the Hamiltonian to the two-dimensional antiferromag- 
netic Ising with coupling constant only in one direction. 
This can be done by the transformations: SVn = 2p, )n — 1 
and pi,„ = l/2(S , i i „ + 1). The Hamiltonian then maps 
onto: 

V \ - 

H = — Si,nSi,n+l (39) 

which has an antiferromagnetic ground state in rows ver- 
tical to chains and therefore the same degeneracy as 
shown above. 

If we turn on the kinetic energy part we observe that 
to first order in perturbation theory, we cannot have non- 
zero matrix elements among the degenerate ground state 

.2 

manifold; then second order processes of order y must 
be taken into account. Therefore the degeneracy is lifted 
by any infinitesimal amount of kinetic energy (tunneling 
strength) which, to second order in perturbation theory, 
selects the checkerboard configuration. It is a version of 
the "order by disorder" phenomenon; order is developed 
by increasing a quantity associated to disorder (kinetic 
energy). To second order in the tunneling strength the 
energy gain is: 



t± t t 

~^ rC i,n C i-n+i C i,n+l C i,n) (41) 

and it favors the checkerboard pattern. This is precisely 
the DW (2kp, tt) that was found using RPA calcula- 
tions. In the pseudospin language, the Ising Hamiltonian 
acquires an additional term which takes into account the 
second order tunneling: 

H = EKj + y)Si,„Si,„+i + (^)s hn s. l+hn }. (42) 

At less than half filling (the physics at greater than 
half filling will be identical through a particle-hole trans- 
formation), the ground state degeneracy in the absence 
of kinetic energy is even worse as there exist extra un- 
occupied sites which could occur anywhere within the 
previous manifold of states. However, now first order 
processes in hopping are allowed. Including the hopping 
along the chain to first order, one can think of the ground 
state as being stacks of particles as delocalized as possi- 
ble, but avoiding contact with other particles above or 
below. For example, at 1/4 filling, each particle will be 
delocalized over two sites in the chains, these stacks then 
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forming an identical problem to the single-site columns in 
the half filled case. Of course at incommensurate fillings, 
the particles will be delocalized over a fractional num- 
ber of sites - but the ultimate physics will remain the 
same. As before, finally taking second-order processes 
into account removes the mass-degeneracy and isolates 
the 'crystalline' phase as the unique ground state. The 
precise determination of the possible favored ground state 
configurations beyond this intuitive argument requires 
more detailed computational work which is beyond the 
scope of the present work and the illustrated physics. 



While such models can still be bosonized in the pres- 
ence of inter-chain backscattering (see e.g. Ref. |55[ ). 
the resulting bosonic Hamiltonian is non-linear and not 
easy to deal with. We therefore elect to remain in the 
fcrmionic representation, where the one-loop RG equa- 
tions can be easily derived: 



dg%, 



dl 

d 92, 



-^gi,nm (<?2, nm ~ .92, nn) ~ ffl.mffl, 



dl 



= -5l,nm- ( 46 ) 



C. The limit of zero interchain hopping 

In the model Eq. Q] for the specific case t± = 0, we 
can study the model in a completely different way, as 
the non-interacting part is now strictly one-dimensional. 
This allows us to single out the most important compo- 
nents of the interaction and look at their flow under the 
renormalization group (RG) rescaling of the system^ 

On each of the chains, one can make the chiral decom- 
position: 



- lkFX L n {i 



a ik F x 



R n {x) 



(43) 



where x — ai, a being the lattice spacing, n indicates 
the chain number and L and R are slowly varying fields 
near the left and right Fermi points respectively. By lin- 
earizing the spectrum around the two Fermi points, the 
kinetic part of the Hamiltonian becomes 

Ho = -ivF^2 / dx ( R i( x )dxR n (x) - L\ l {x)d x L n {x)) . 

where vf is the (bare) Fermi velocity, and the interaction 
term becomes 

nm 

where g2,nm is the forward scattering, and gi, nm is the 
back scattering between particles on chains n and m. 
In terms of the bare parameters of the Hamiltonian, 
<?i,nm — 52, nm = V for nm nearest neighbors and zero 
otherwise - however under RG, effective further neighbor 
scattering terms are generated, so we must retain them in 
the Hamiltonian. We note that there is also inter-chain 
interactions involving only right (or left) moving particles 
- however these are known not to affect the RG flow^ so 
we will not consider them. 

If only the forward scattering (g^) terms were present, 
this model can be exactly solved by the Bosonization 
technique.— There is a transition to a smectic phased at 
large V ~ vp. As this transition will only be visible in a 
slightly perturbed version of the model described by Eq.Q] 
with longer range interactions 39 which we have previously 
discussed, we will concentrate here on the effect of the g\ 
backscattering terms. 



where / = —In A, with A being the ultra-violet cut-off 
of the theory. These are very similar (although not 
identical) to the related problem of coupled spin-full 
chains 

Following Ref. 571, we can now see that while the bare 
couplings may be nearest neighbors only, the RG due to 
the sum in Eq. 1461 generates couplings between all pairs of 
chains. Furthermore, due to this sum, the backscattering 
amplitudes renormalize much quicker than their forward 
scattering counterparts^, so to leading order one can 
neglect the renormalization of gi , and deal with the far 
simpler equation: 



dgi 



dl 



i9U 



(47) 



We note that effects of 32 were dealt with in Ref. [58| 
which concludes that the only interesting effects of 32 
come when one has a bare intra-chain backscattering am- 
plitude, something we neither have by construction of our 
interaction, nor can have due to the spinless nature of the 
model. We therefore will deal only with Eq. |47l which 
are easily separated by a Fourier transform 



9i 



i(n-m)q ± 
yi.nm 



(48) 



which gives us 



9gi(q±) 
dl 



-g\{g± 



(49) 



Hence any parts of gi(q±) which begin negative renor- 
malize to strong coupling. Now, the bare backscattering 
amplitude from nearest neighbors is 



gr e (q±) = 2Fcos(< 7x ), 



(50) 



so under RG, the component of the interaction which 
hits strong coupling first (after which the RG must be 
stopped) is the q± = n part of the backscattering ampli- 
tude, or in other words, the (2fc^,7r) component of the 
interaction. This signifies the formation of a (2kp,Tr) 
DW 5 iS as we expected. 
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v/t 



FIG. 8: [Color online] Plot of the density-wave and meta- 
nematic critical surfaces for fi/tn = —1.5. The plots are qual- 
itatively similar for other values of \i. 



IV. FULL PHASE DIAGRAM AND EFFECT OF 
TRAPPING POTENTIAL 

A. Full phase diagram 

We now can draw a complete phase diagram of the 
model, including both possible instabilities - which is 
plotted in Figure [8] for a typical value of /z/tii = —1.5. 
The plots are qualitatively similar for other values of 
the chemical potential. From the plot, it may seem that 
the instability to a crystalline pattern dwarfs the meta- 
nematic instability, however it is also clear that there will 
always be a region in phase space where the meta-nematic 
transition can be observed. Furthermore, we should add 
that while in our simple model, the density-wave is usu- 
ally the dominant instability, the two transitions come 
from fundamentally different physics: the meta-nematic 
transition from the van-Hove singularity in the density 
of states, and the density-wave transition from a near- 
nested Fermi surface. Therefore, destroying the nesting 
by for example adding a next-nearest-neighbor hopping 
term; or frustrating the interaction will both lead to a 
suppressed tendency to enter the DW phase whilst barely 
affecting the meta-nematic one. However, the best way 
to enhance the MN phase relative to the DW is open for 
further investigation. 

As noted in the Introduction the separation in phase 
space between the MN and DW phase transitions is one 
of the salient features of this particular model, with the 
phase diagram in Fig. [5] showing clearly the existence of 
such regions of phase space where either one or the other 
instability takes place. On the other hand, in those other 
regions where the two different instabilities become too 
close together, our calculational techniques break down, 
as the mutual competition between each of them must 
be carefully taken into account. We discuss this point 
further in Section V below. For the time being we simply 
must bear in mind that the plotted phase diagram is less 



reliable around the latter regions. 



B. Density profiles in a trapping potential 

Having discussed the properties of transitions in a gen- 
eral context, we now turn specifically to one potential re- 
alization of the model: dipolar fermions in an anisotropic 
optical lattice.— In a cold atoms experiment, in addition 
to the optical lattice, there is also a harmonic trapping 
potential: 

V(r)=ar 2 , (51) 

where a is some constant dependent on the experimental 
setup. Within the local density approximation (LDA), 
this can be theoretically approximated as an inhomoge- 
nous chemical potential 

Meff(r) =n-V(r) = [i-ar 2 . (52) 

From the previous formulas, we can therefore plot the 
density profile, and profile of renormalized hopping) 
within the trap, which is shown in Fig. |H1 We note that 
in this Figure, we show only the metanematic transition 
and not any possible DW phases: in fact for the rather 
strong interaction strengths shown in the figure the DW 
phase will have a dominant presence in the phase dia- 
gram. However on reducing V, the DW will go away, 
leaving only the metanematic jump shown. This jump 
will be present for any interaction strength V, although 
for small V , the size of the jump becomes exponentially 
small, which is why it is artificially enhanced in the fig- 
ure, by plotting the profile for large values of V (see also 
the discussion in the previous section about the relation 
between the DW and MN transitions). 

The interesting feature of these steps in the density 
profile, is that they signal the separation of two compress- 
ible phases,— and do not represent the more common 
compressible/incompressible boundary. The only differ- 
ence between the phases on cither side of the jump is the 
topology of the Fermi surface: on one side it is closed, 
while on the other it is open. This difference is, how- 
ever, possible to measure in optical lattice experiments 
by direct imaging of the Fermi surfacei 61 ' 62 



V. DISCUSSION 

We have examined in detail two different phase tran- 
sitions associated with quantum liquid crystal phases of 
matter, which appear in the same physically realizable 
model. The first, the meta-nematic transition which 
is non-symmetry breaking but has strong ties to both 
the Pomeranchuk and the Lifshitz transitions while the 
second crystallization transition is a standard symmetry 
breaking transition. We have shown that both transi- 
tions may be found in a simple anisotropic lattice model 
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FIG. 9: [Color online] The local density as a function of po- 
sition of our model within a harmonic trap (treated within 
the LDA). The parameters are adjusted so that the density 
in the middle of the trap is exactly half filled, and the bare 
transverse hopping — 0.1. We note that the variation 
of the effective chemical potential across the trap makes the 
meta-nematic transition visible as a jump in the density, with 
no special tuning of parameters. Inset: the variation of the 
'order parameter' tj_ aginst position in the trap. 



of spinless fermions, which may be realized with ultra 
cold dipolar atoms in an optical lattice. 

The Pomeranchuk mechanism of spontaneous de- 
formation of the Fermi surface is a straight-forward 
generalisation of the Stoncr mechanism of itinerant 
ferromagnetismi (very much in the same way as un- 
conventional pairing generalises s-wave pairing in the 
BCS theory of superconductors). This begs the question 
why, unlike itinerant ferromagnetism, which is ubiqui- 
tous, there is so little evidence of Pomeranchuk and other 
spontaneous deformations of the Fermi surface, such as 
the meta-nematic transition discussed here, in real sys- 
tems (in contrast, many unconventional superconductors 
exist in addition to those conforming to the original form 
of BCS theory). One insight is offered by the mean field 
analysis of central interactions leading to Pomeranchuk 
instabilities^ It reveals that unlike the Stoner instabil- 
ity, which may be induced by a contact interaction, its 
higher-angular momentum generalisations, namely the 
Pomeranchuk instabilities, require a repulsive interaction 
potential V(r) with non-trivial features at a characteris- 
tic distance r > r s , where r s is the average distance 
between particles. Under such conditions other insta- 
bilities, involving the breaking of translation symmetry, 
are also expected. In particular there is a competition 
between Fermi surface deformations and more conven- 
tional, density-wave instabilities leading to smectic order 
(stripes) or crystallisation. A similar situation is found 
in mean field theories of Fermi surface deformations on 
the square lattice : 63 ' 64 it is found that repulsion between 
fermions sitting on nearest-neighbour sites can lead to 
Pomeranchuk deformations of the Fermi surface but that 



this tendency is strongest near half filling where other, 
real-space instabilities are also enhanced. The question 
of which type of order prevails in this case is hard to 
answer. In general, it will be model-specific and require 
a quantitative asessment of the role of fluctuations go- 
ing beyond mean field theories. It is worth mentioning 
in passing, that recently, there was a rigorous study of 
the effects of van Hove points on the properties of in- 
teracting fermionic systems, using functional renormal- 
ization group (FRG) techniques. 65 There it was shown 
that the logarithmic singularity in the density of states 
caused by the zero of the bare Fermi velocity gives rise 
to new marginally relevant terms in the renormalization 
group equations and leads to singularities in the analytic 
properties of electronic self-energies. 

In the present work we have considered an anisotropic 
square lattice with anisotropic interactions and we have 
found that in this case the situation is much more clear- 
cut. As discussed above, the density wave and meta- 
nematic instabilities are driven by two separate aspects 
of the band structure, namely nesting and the van Hove 
singularity in the density of states, respectively As is well 
known, these two features coincide for a square lattice, 
but as soon as t± ^ tn they diverge: while perfect nesting 
still occurs at exactly half filling, // = 0, the singularity 
splits into two separate features at fi = ±2|tii —t±\. As a 
result of this, for the anisotropic model considered here 
there is a broad region of the phase diagram where meta- 
nematic order clearly prevails over the tendency to form 
a density wave. 

This separation of features also partially justifies the 
use of the Hartree-Fock approximation - far enough away 
from the density- wave instability, the susceptibility of the 
model remains small, so the interaction is only weakly 
renormalized by higher order terms, leading to only small 
corrections to the Hartree-Fock results. As the two dif- 
ferent instabilities come closer together in phase space 
(as would always happen near the Pomeranchuk insta- 
bility in lattice systems), the interplay between the FS- 
deforming transition and the crystalline transition must 
also be considered, for example using the fluctuation- 
exchange approximation (FLEX). 66 A closely related 
competition between spin density wave formation and su- 
perconductivity has recently been studied in the context 
of the pnictide superconductors^ Such a study however 
is beyond the scope of the present work. 

We also note that even far from the crystalline transi- 
tion, the Hartree-Fock approximation neglects all fluctu- 
ations, and it is natural to ask what role such fluctuations 
may play in the meta-nematic transition. We expect that 
as this transition is driven by a van-Hove singularity in 
the density of states, fluctuations which can be described 
by Fermi-Liquid physics will not qualitatively change the 
picture. The reason is that so long as the energy distri- 
bution function at T — has a jump at the Fermi energy, 
the same logarithmic divergences as we have considered 
will appear in the perturbation series, which will lead to 
the same physics, and the effect of the fluctuations will 
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be limited to the renormalization of some parameters. 

Interestingly, this separation of features due to 
anisotropy is somewhat reminiscent of the cuprate 
high-temperature superconductors, where the exis- 
tence of stripe order is well-established in tetragonal 
Lai.4gNdo.4Sro.i2Cu04 and Lai.875Bao.i2sCu04 [Ref. lo8l 
and references therein] while evidence is accumulating of 
a nematic state in YBa2Cu306 +a; i 69 i 70 The CuO chains 
present in the latter material constitute a very strong or- 
thorhombic perturbation. One way of interpreting the 
date is that the resulting breaking of the D^h symmetry 
of the crystal lattice drives away a density wave instabil- 
ity that competes with the nematic order, allowing the 
latter to prevail in the case of YBCO. 

We finally note that in the general case when interac- 
tions between next nearest neighbours on different chains 
become important the density wave instability described 
here may take place with wave vector Q = (0, n) (lead- 
ing to a stripe phase), instead of (2kp, ir) (checkerboard). 
The stripe instability becomes preferable due to a frus- 
tration of the nearest neighbour interaction^ 9 Similarly 
in some cuprates a checkerboard pattern is observed in- 
stead of stripes or nematic order [Ref. [7l|, and references 
therein] . 

In summary, we have studied in detail a toy model 
which shows a very rich phase diagram and is a play- 
ground for many phases which are sought after by the 
community of correlated fermions. We have demon- 
strated within this model the fundamentally different 
properties of the Fermi surface that lead to crystalliza- 
tion and a Fermi-surface-modifying transition. We must 
add that while our analysis shows the presence of a meta- 
nematic transition, it is exponentially small in the inter- 
action strength, and therefore presents a formidable task 
at the moment to measure it experimentally. At larger 
interactions of course, the perturbation theory needs fur- 
ther supplementary work by treating both the meta- 
nematic transition and the density wave instability on 
the same footing with e.g. FLEX, or with strong cou- 
pling techniques such as FRG. We hope that this work 
will stimulate more studies in this direction, both exper- 
imentally and theoretically. 
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Appendix A: Derivation of the free energy 
expansion near the quantum critical endpoint 

In this Appendix, we derive the expansion Eg. 1191 from 
the free energy expression Eq. [T5] at T = 0. While the 
energy itself is a continuous function, the presence of van- 
Hove singularities in the density of states leads to loga- 
rithmic divergences in derivatives. The goal is as follows, 
following the spirit of Landau: to expand the Free en- 
ergy near the van-Hove energies to obtain the leading 
order logarithmic terms, and then by using phenomeno- 
logical arguments in performing the rest of the calcula- 
tions and adding the different terms while keeping the 
correct physics, we obtain the expansion Eq. [T5] in the 
main text. 

The kinetic energy can be written as (we measure en- 
ergies in units of t\\) 

(Fo) = J ^£ (k)n(k) = 9 2(x - x Q )f (Al) 
where 

9 = J (g^(kWk) 

/ = / 70jcos^n(k). (A2) 

Some elementary manipulations transform the potential 
energy into the form 

= -Vf. (A3) 
We now define the normal density of states 

J .0^6 (t- em, (A4) 

and a weighted density of states 

*(0=2tt J J^cosM(£-£*(k)). (A5) 
This allows us to write our integrals as 




where —A is the bottom of the band. As we will be 
interested in the derivative of the energy with respect to 
small changes in the distribution function, the values of 
such energy integrals at their lower limit is unimportant, 
and will be ignored in the next. 

The density of states has a logarithmic divergence 
at the van-Hove energies £ v h , which is a result of the inte- 
gration around the saddle points at (±7r, 0) and (0, ±7r). 
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While the energy of these two points coincides in the 
case t± — tii, for the anisotropic model we study, they 
occur at two different energies. For the case t± < ti and 
less than half filling (which is what we deal with in the 
present manuscript), the important points are (0, ±7r). 
This gives 

^)oc-ln(C-^H) (A7) 

and similarly as cosk y = — 1 at these points, 

HO « ln(£ - Cvh). (A8) 

Taking into account that the van-Hove energy is pro- 
portional to the parameter x, it transforms the leading 
logarithmic terms into 



g oc x 2 In | x \ + regular terms 
/ oc x In | x | + regular terms. 



At this stage, while the regular terms are sub-leading and 
thus unimportant for small enough x, a little care must 
be taken with them in order to obtain the correct physical 
behavior. In particular, the kinetic energy g — 2(x — Xo)f 
whilst showing the above non-analytic behavior as x — > 
must also have a minimum at x = xq. The simplest 
addition of an x 2 term ensures this, giving 



(Fo) 



\n\x\-- 



2x [x — Xo] In \x\ 



(A10) 



(A9) 



As the potential energy is proportional to f 2 , one must 
here also consider the constant term in / ~ c + a; In \ x\, 
which gives the expansion Eq. [TO] in the main text . 



1 I. Ia .Pomeranchuk, Sov. Phys. JETP 35, 524 (1958). 

2 H. Yamase and H. Kohno, J. Phys. Soc. Japan 69, 2151- 
2157 (2000). 

3 C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162- 
5165 (2000). 

4 B. Valenzuela and M. A. H. Vozmediano, Phys. Rev. B 63, 
153103 (2001). 

5 V. Oganesyan, S. A. Kivelson and E. Fradkin, Phys. Rev. 
B 64, 195109 (2001). 

6 R. Settai et al, J. Phys. Soc. Jpn. 74, 3016 (2005). 

7 H. Sugawara et al, Phys. Rev. B 66, 134411 (2002). 

8 S. Paschen et al, Nature 432, 881 (2004). 

9 K. B. Cooper, M. P. Lilly, J. P. Eisenstein, L. N. Pfeiffer 
and K. W. West Phys. Rev. B 65, 241313(R) (2002). 

10 S. A. Grigera et al, Science 306, 1154 (2004). 

11 J.-F. Mercure et al, Phys. Rev. Lett. 103, 176401 (2009). 

12 J. Quintanilla and A. J. Schofield, Phys. Rev. B 74 115126 

(2006) . 

13 C. M. Varma and L. Zhu, Phys. Rev. Lett. 96 36405 (2006). 

14 C. M. Varma and L. Zhu, Phys. Rev. Lett. 98 177004 

(2007) . 

15 Q. M. Doan and E. Manousakis, Phys. Rev. B 75 195433 
(2007). 

16 H. Yamase, Phys. Rev. Lett. 102, 116404 (2009); 

17 H. Yamase, Phys. Rev. B 80, 115102 (2009). 

18 Wei-Cheng Lee and Congjun Wu, Phys. Rev. B 80, 104438 
(2009). 

19 M. Zacharias, P. Wolfle, and M. Garst, Phys. Rev. B 80, 
165116 (2009). 

20 C. M. Puetter, J. G. Rau, and Hae- Young Kee, Phys. Rev. 
B 81, 081105(R) (2010). 

21 N. Shannon, T. Momoi and P. Sindzingre, Phys. Rev. Lett. 
96, 027213 (2006). 

22 S. A. Kivelson, E. Fradkin and V. J. Emery, Nature 393, 
550 (1998). 

23 E. Fradkin, S. A. Kivelson and V. Oganesyan, Science 315, 
196 (2007). 

24 C.H. Chen and S.-W. Cheong, Phys. Rev. Lett. 76, 4042 
(1996). 



25 J. M. Tranquada, B. J. Sternlieb, J D. Axe, Y. Nakamuara 
and S. Uchida, Nature 375, 561 (1995). 

26 D. G. Clarke, S. P. Strong and P. W. Anderson, Phys. Rev. 
Lett. 72, 3218 (1994). 

27 E. Arrigoni, Phys. Rev. Lett. 83, 128 (1999). 

28 E. Arrigoni, Phys. Rev. B 61, 7909 (2000). 

29 S. Ledowski and P. Kopietz, Phys. Rev. B 76, 121403 
(2007). 

30 J. J. Betouras and J. T. Chalker, Phys. Rev. B 62, 10931 
(2000). 

31 V. L. Campo, K. Capelle, J. Quintanilla and C. Hooley, 
Phys. Rev. Lett. 99, 240403 (2007). 

32 Tin-Lun Ho and Qi Zhou, Nat. Phys. 6, 131-134 (2010). 

33 M. A. Baranov, L. Dobrek, and M. Lewenstein, Phys. Rev. 
Lett. 92, 250403 (2004). 

34 H. P. Biichler, E. Dernier, M. Lukin, A. Micheli, N. 
Prokof'ev, G. Pupillo and P. Zoller, Phys. Rev. Lett. 98, 
060404 (2007). 

35 A. Micheli, G. Pupillo, H. P. Biichler and P. Zoller, Phys. 
Rev. A 76, 043604 (2007). 

36 K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, B. 
Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, D. 
S. Jin, and J. Ye, Science 322, 231 (2008). 

37 Mingwu Lu, Seo Ho Youn, and Benjamin L. Lev, Phys. 
Rev. Lett. 104, 063001 (2010). 

38 T. Miyakawa, T. Sogo, and H. Pu, Phys. Rev. A 77, 
061603(R) (2008). 

39 J. Quintanilla, S. T. Carr, J. J. Betouras, Phys. Rev. A 
79, 031601(R) (2009). 

40 S. T. Carr, J. Quintanilla, J. J. Betouras, Int. J. Mod. 
Phys. B 23, 4074 (2009). 

41 B. M. Fregoso, K. Sun, E. Fradkin and B. L. Lev, New J. 
of Physics, 11, 103003 (2009). 

42 N. R. Cooper, and G. V. Shlyapnikov, Phys. Rev. Lett. 
103, 155302 (2009). 

43 Y.-P. Huang and D.-W. Wang, Phys. Rev. A 80, 053610 
(2009). 

44 C.Zhao, L.Jiang , X.Liu, W.M.Liu, X.Zou and H.Pu, 
larXiv:0910.4"775 (2009). 



14 



45 C.Lin, E.Zhao and W.V.Liu, Phys. Rev. B 81, 045115 
(2010). 

46 I. M. Lifshitz, Sov. Phys. JETP 11, 1130 (1960). 

47 R. S. Perry et al, Phys. Rev. Lett. 86, 2661 (2001). 

48 S. A. Grigera et al, Science 294, 329 (2001). 

49 Y. Yamaji, T. Misawa and M. Imada, J. Phys. Soc. Jpn. 



50 



75, 094719 (2006). 

Y. Okamoto, A. Nishio, and Z. Hiroi, Phys. Rev. B 81, 
121102(R) (2010). 

51 D. M. Edwards, Physica B 378-380, 133-134 (2006) 

52 y Yamaji, T. Misawa and M. Imada, J. Phys. Soc. Jpn. 

76, 063702 (2007). 

63 I.E.Dzyaloshinskii and A.I.Larkin, Sov. Phys. JETP, 38, 
202 (1974) 

54 J. Solyom, Advances in Physics 28, 201 (1978). 

55 S. T. Carr, and A. M. Tsvelik, Phys. Rev. B 65, 195121 
(2002). 

56 L. P. Gorkov and I. E. Dzyaloshinskii, Zh. Eksp. Teor. Fiz. 
67, 397 (1974). 

57 L. Mihaly and J. Solyom, J. Low Temp. Phys. 24, 579 
(1976). 

58 N. Menyhard, Solid State Communications 21, 495 (1977). 



59 J. Solyom, J. Low. Temp. Phys. 12, 547 (1973). 

60 Yong-il Shin et al, Nature 451, 689-693 (2008). 

61 I. Bloch, Nature Physics, 1, 23 (2005). 

62 M. Kohl, H. Moritz, T. Stoferle, K. Giinter and T. 
Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 

63 B. Valenzuela, and M. A. H. Vozmediano, Phys. Rev. B 
63, 153103 (2001). 

64 J. Quintanilla, C. Hooley, B.J. Powell, A.J. Schofield, and 
M. Haque. Physica B 403, 1279 (2008). 

65 J. Feldman, M. Salmhofer, Rev. Math. Phys. 20, 275 
(2008). 

66 N. E. Bickers, D. J. Scalapino and S. R. White, Phys, Rev. 
Lett. 62, 961 (1989). 

67 A. V. Chubukov, D. V. Efremov and I. Eremin, Phys. Rev. 
B 78, 134512 (2008). 

68 P. Abbamonte, A. Rusydi, S. Smadici, G. D. Gu, G. A. 
Sawatzky, and D. L. Feng, Nature Physics 1 155 (2005). 

69 Y. Ando, K. Segawa, S. Komiya, A. N. Lavrov, Phys. Rev. 
Lett. 88, 137005 (2002). 

70 V. Hinkov et al, Science 319, 597 (2008). 

71 W. D. Wise et al, Nature Physics 4, 696 (2008). 



